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The Soft Function for Exclusive iV-Jet Production at Hadron Colliders 
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The iV-jettiness event shape divides phase space into N + 2 regions, each containing one jet or 
beam. Using a geometric measure these regions correspond to jets with circular boundaries. We give 
a factorization theorem for the cross section fully differential in the mass of each jet, and compute 
the corresponding soft function at next-to-leading order (NLO). The ultraviolet divergences are 
analytically extracted by exploiting hemispheres for interactions between each pair of hard partons, 
leaving only convergent integrals that are sensitive to the precise boundaries. This hemisphere 
decomposition can also be applied to other iV-jet soft functions, including other observables. For 
iV-jettiness, the final result for the soft function involves stable one-dimensional numerical integrals, 
and all ingredients are now available to extend NLO cross sections to resummed predictions at 
next-to-next-to-leading logarithmic order. 



I. INTRODUCTION 

The measurement of exclusive jet cross sections, where 
one identifies a certain number of signal jets but ve- 
toes additional jets, is an important aspect of Higgs and 
new-physics searches at the LHC and Tevatron. In such 
searches, the experiments often analyze the data sepa- 
rated into bins of different numbers of jets. This is done 
because the relative contributions of various signal and 
background channels often vary with the number of hard 
jets in the event. Hence, it is important to optimize the 
analyses for each jet bin. An important example is the 
Higgs search at the Tevatron [lj , which analyzes the data 
separately for H + jets, H + 1 jet, and H + 2 or more 
jets. 

Reliable theoretical calculations of exclusive jet cross 
sections are of course essential. The complication com- 
pared to the calculation of an inclusive iV-jet cross sec- 
tion, where one sums over any additional jets, comes from 
the fact that the veto on additional jets imposes a restric- 
tion on the energetic initial- and final-state radiation off 
the primary hard partons, as well as the overall soft ra- 
diation in the event. This restriction on additional emis- 
sions leads to the appearance of large Sudakov double 
logarithms in the perturbation theory. This is a well- 
known phenomenon and is due to an incomplete cancel- 
lation of infrared contributions between virtual correc- 
tions and restricted real radiation. For this reason, the 
calculation of exclusive jet cross sections is traditionally 
carried out with parton-shower Monte Carlo programs, 
where the parton shower allows one to resum the most 
singular leading double logarithms. 

An alternative analytic approach to calculate exclusive 
jet cross sections is possible using factorization and the 
methods of soft-collinear effective theory (SCET) 
SCET is designed to study processes with a specific num- 
ber of hard jets. It allows one to factorize the iV-jet cross 
section into individually calculable pieces and resum the 
large logarithmic contributions to obtain a convergent 



perturbative series. The advantage of this approach is 
that the resummation can be carried out to much higher 
orders than is possible with parton showers. In addition, 
it is much easier than in a parton-shower program to in- 
clude higher-order virtual corrections, and to correctly 
reproduce the inclusive cross section in the limit when 
the jet- veto cut is eliminated. 

Schematically, the cross section for pp — > N jets (plus 
some nonhadronic final state which we suppress for now) 
can be factorized as 
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This formula directly applies to observables that im- 
plement a veto on additional jets which restricts the 
phase space to the exclusive iV-jet region (assuming that 
Glauber effects cancel as they do in Drell-Yan Q). The 
hard function Hn encodes hard virtual corrections to 
the underlying partonic 2 — > N process, the beam func- 
tions B a i, contain the parton distributions and pertur- 
bative collinear initial-state radiation from the colliding 
hard partons, and the jet functions Ji describe energetic 
collinear final-state radiation from the primary N hard 
partons produced in the collision. The soft function Sn 
describes the soft radiation in the event that couples to 
the in- and outgoing hard partons. Since the collinear 
and soft radiation are not separately observable, the soft 
function is convolved with the beam and jet functions. 
The veto on additional jets restricts the collinear initial- 
state radiation, the final-state radiation, and the soft 
radiation, which means the precise definition of the re- 
quired beam, jet, and soft functions depends on the veto 
variable. 

For the case of an exclusive 0-jet cross section, inclu- 
sive beam and jet functions can be obtained by using 
a simple event-shape variable called beam thrust Q to 
veto central jets. This 0-jet cross section has been stud- 
ied for Drell-Yan production in Ref. [l(J and for Higgs 
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production in Ref. jllj. The latter is for example rele- 
vant for the H — > WW* search channel, where a jet veto 
is needed to remove the large tt — > WWbb background. 
The use of an event shape for the jet veto makes a re- 
summation of large logarithms to next-to-next-to- leading 
logarithmic (NNLL) order possible. 

The generalization of beam thrust to processes with 
N jets is iV-jettiness, Tn, introduced in Ref. [T^- It 
is designed such that in the limit Tn — > the final state 
consists of N narrow jets plus two narrow initial-state ra- 
diation jets along the beam axis (for hadron collisions). 
Since it does not restrict the collinear radiation inside a 
jet, the beam and jet functions appearing in Eq. ^ arc 
again the inclusive beam and jet functions (which are 
known to one [l(| El El El and two loops EH, re- 
spectively). Furthermore, since 7V-jettiness itself covers 
all of phase space, no additional restriction on the ra- 
diation outside of jets or beams is needed. In contrast, 
hadron-collider event shapes constructed from transverse 
momenta only, such as transverse thrust, in general re- 
quire the addition of exponentially suppressed forward 
terms to suppress the contributions from large rapidi- 
ties [13, El- 

Factorization for iV-jettiness can be contrasted with 
factorization for jet algorithms. Here, the perturbativc 
corrections are complicated by: the presence of non- 
global logarithms |19l - t22l | , the potential for soft radiation 
to be strongly influenced by the number of energetic par- 
tons in the jets, and by cuts on soft radiation that intro- 
duce additional soft scales that must be handled within 
factorization 0, [24| . Jet functions for jet algorithms in 
e + e~ — > jets have been calculated at next-to-leading or- 
der (NLO) in Refs. (HH^. In Ref. [H the soft function 
for e + e~ —> jets was calculated at NLO, where a cut on 
the total energy outside the jets was used as the jet veto. 
Using A^-jettiness avoids these issues that complicate the 
structure of perturbation theory. 

The iV-jettiness event shape assigns all particles to one 
of 7V+2 regions, corresponding to the N jets and 2 beams. 
Therefore Tn acts much like a jet algorithm, and we can 
consider distinct measurements on each of these "jets" . 
The simplest example is T N , the A^-jettiness contribution 
from each region i, where Tn — Tiv- A measurement 
of T N is essentially the same as measuring the transverse 
mass or invariant mass of this jet. This correspondence 
will be made precise in the next section. We will also 
briefly explore the shape of the jet regions obtained using 
iV-jettiness with different measures. A geometric mea- 
sure gives jets with circular boundaries, putting them in 
the class that are typically preferred experimentally. 

For an iV-jettiness cross section calculation using 
Eq. (TT]), the only missing ingredient for an evaluation 
of generic processes at NNLL is the one-loop iV-jettiness 
soft function, Sn, which we compute in detail in this 
paper. (As mentioned above, the beam and jet func- 
tions are known. The hard function in Eq. ([1} can be 
obtained from the corresponding QCD fixed-order cal- 
culation, many of which are now known to NLO.) Gen- 



eral features of A^-jettiness and its jet regions are ex- 
plored in Sec. HU Results are given for the fully differ- 
ential T N factorization theorem, and for renormalization 
group consistency equations for the A^-jettiness soft func- 
tion. Section IIIII contains details of the NLO calculation 
of Sn, including developing a simple method that uses 
hemispheres for each pair of hard partons to extract UV 
divergences and the corresponding induced logarithmic 
terms. The remaining 0(a s ) terms are then given by fi- 
nite integrals that do not involve the UV regulator, and 
we will refer to these as the non-hemisphere contribu- 
tions. This hemisphere decomposition is not specific to 
the iV-jettiness observable, and we show how it can be 
applied in general. For the iV-jettiness soft function we 
reduce the non-hemisphere contributions to well-behaved 
one-dimensional numerical integrals (some details are rel- 
egated to appendices). Section HVl contains conclusions. 

Although it is not directly related to our investigations 
here, it is worth mentioning that iV-jettiness is useful for 
exploring jet substructure j2(|[23]. This is done with N- 
subjettiness, which restricts the definition of the event 
shape to particles and reference momenta inside a jet. 
There are interesting correspondences between applica- 
tions of iV-jettiness and ./V-subjettiness. In particular 
one could study the mass spectrum of subjets with T N , 
following a similar procedure that we advocate here for 
jets. 



II. SETUP OF THE CALCULATION 
A. AT-Jettiness Definition and Regions 

A^-jettiness is defined as [lH 

where i runs over a, b for the two beams and 1, . . . , N for 
the final-state jets. 1 For e + e~ collisions, the terms for 
the beams are absent and we continue to let A~ refer to 
the number of jets. The complexity of the calculation for 
the e + e _ (N + 2)-jettiness is equivalent to A^-jettiness 
for pp collisions. In Eq. © the qi are massless reference 
momenta for the jets and beams, and the Qi are normal- 
ization factors. For each jet we can take 

gf = uji (1, Hi) , (3) 

where u>i is the jet energy, and Hi is the jet direction. 
The uji and Hi can be predetermined with a suitable 
jet algorithm, and the choice of algorithm o nly gives 
power-suppressed effects, as explained in Ref [12]. For 



1 Here we use a dimension-one Tn as in Ref. [Tlj as opposed to 
the dimensionless rjv of Ref. \\'M . 
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FIG. 1: Jet and beam reference momenta for 1-jettiness (left), 2-jettiness (middle) and e + e 3-jettiness (right). In the middle 
plot the jets and beams do not necessarily lie in a plane. 
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FIG. 2: The jet and beam regions for the same two jets using 2-jettiness. On the left we use the invariant-mass measure 
Qi = Q. On the right we use the geometric measure with Qi — \qir \ for the jets and Q a ,b = x a:b E cln for the beams. 



the beams we have 

Qa = \ x a E cm (l, z), qg = ix fc E cm (l, -z) , (4) 

where E cm is the center-of-mass energy, z points along 
the beam axis, and x a-b are the light-cone momentum 
fractions of the colliding hard partons. The latter are 
defined as 

x a E cm = Q e Y , x b E cm = Q e~ Y , (5) 

where Q 2 and Y are the total invariant mass-squared and 
rapidity of the hard interaction. They are determined 
from the observed final state by 

Q 2 = x a x b E 2 m = (qi H h q N + q) 2 , 

2r = ln^=ln ( ^~f ) ; (gl + -- + gJV + ; ) . (6) 
x b (1, z) ■ (qi H h qN + q) 

Here q^ denotes the total momentum of the non-hadronic 
final state if one is present. 

The choice of the q£ is illustrated in Fig. [T] for 1- 
jettiness (left panel), 2-jettiness (middle panel), and 
e+e~ 3-jettiness (right panel). For the first two cases 
q^ is given by the momentum of the W/Z. In SCET 
the qtf's become the large label momenta on the collinear 
fields, which can be thought of as the momenta of the 
partons in the hard interaction. The minimum in Eq. @ 



divides the total phase space into N + 2 regions, one for 
each beam and jet, as indicated by the dashed lines in 
Fig.[U Their union exactly covers all of phase space, and 
the boundary between any two regions is a (part of a) 
cone. 

The Qi in Eq. @ are dimension-one variables that 
characterize the hardness of the jets. Different choices 
for the Qi correspond to choosing different distance mea- 
sures in the minimization in T/v- For example, for fixed 
Qi = Q, the distance measure is just the invariant mass, 
2gj • pk- The resulting jet and beam regions in this case 
are illustrated for 2-jettiness in the left panel of Fig. [21 
Choosing the jet transverse momentum Qi = \qir \ for the 
jets, the measure becomes a geometric measure, which is 
boost-invariant along the beam axis, 

2 1l P , k = \Pkr\ (2coshA?7 ife - 2cosA0 lfe ) 

\QiT\ 

w \p kT \ [(AtmO 2 + (A<^) 2 ] . (7) 

Here, Arjik = r\i — r)k, &4>ik — 4>%~ 4>k are the differences 
in (pseudo) rapidity and azimuthal angle between the di- 
rection of jet i and particle k. The second line is valid 
in the limit of small A77 and Acf>. Equation results 
in circular boundaries for the jet regions, as illustrated 
in the right panel of Fig. [5] In this case only the Hi 
part of gf enters, and the Hi could be obtained by the 
choice which minimizes T/v, thus making ./V-jettiness a 
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true event shape that does not depend on any auxiliary 
input from a jet algorithm. The jet energy is then simply 
given by summing over the particles in each jet region as 
determined by Tn- 

For the beams we have 



2gq ■ Pk = Q e r- 

Qa Qa 



(8) 



with Y — r)k — > —Y + rjk for a — ¥ b. Here two potential 
choices for Q a ^ b are Q a ,b — Q, giving the invariant-mass 



distance measure, or Q a ,b 
gives 

2gq,b ■ Pk _ 
Qa,b 



= Qe 



,±Y _ 



,E rm , which 



\Pkr\e 



(9) 



We will carry out our analysis and one-loop calcula- 
tions keeping Qi arbitrary, enabling various choices to be 
explored using our results. From an experimental point 
of view certain choices will be more advantageous than 
others. For example, the second choice above for Q a ,b is 
useful if the total rapidity cannot be measured because 
there are missing-energy particles in the final state. 

For convenience we define the dimensionless reference 
momenta and their invariant products 



Sij 2f/j * Qj 

We can then rewrite Eq. © as follows, 



(10) 



i i 

P t = J2PkIi e (^-Pk-qi-Pk), 



(11) 



where Pf is the total four- momentum in region i. The 
T N are thus given by the small light-cone component of 
the Pi measured along their respective collinear direc- 
tions qi. In the next section we explore the factorization 



theorem that is fully differential in the T N - The result- 
ing fully differential soft function will be the focus of our 
calculations. 



B. iV-Jettiness Differential in Jet Regions 

The factorization theorem for dtr/dTzv was given in 
Ref. [IH , and is derived in a straightforward manner from 
SCET, see Refs. (with an assumption so far implicit 
in all iV-jet factorization formulae about the cancellation 
of Glauber gluons) . Instead of measuring Tn > the manip- 
ulations leading to the factorization theorem are no more 
difficult when we consider the fully differential cross sec- 
tion, where we measure each individual T N - The value 
of T N determines the transverse mass of region i relative 
to the direction Hi since 



M$ T = P( 



2q 

QiU [l 



-Pix = (fii ■ Pi)(m ■ Pi) 
P[l + 0(\ 2 )] 



o{\ 2 )} 



(12) 



where 



= (1,-nj 



In the last line we 



used fii ■ qi = fii ■ Pi + 0(X 2 ), where A 2 ~ T N /Q and the 
power corrections depend on how the magnitude of % is 
fixed. 

If the label vector Hi is chosen to be aligned with 
the direction of the jet three-momentum Pi such that 
Hi ■ P/\Pi\ ~ 1 + C(A 4 ) then P? ± = + C(A 4 ) and the 
transverse mass is the same as the invariant mass. 



Mf = Pf = Q i T N [l + 0{X 2 )] 



(13) 



Thus the differential T N spectrum corresponds to the 
spectrum in the invariant mass for jet i, where M 2 — > 
for a pencil like jet of massless partons. 

The factorized form for the cross section in the limit 
where all the T N are assumed to be parametrically com- 
parable but small compared to Qi ~ Q is 



da 



AT^dT N ---dTff 



= Jdx a dx b Jd i qd<5> L {q) J d<S> N ({q.j}) M N (<f> N , $l) (27r) 4 <5 4 (<7 a + q b - q% q N ~ q) 

N 

x ^ dt a B Ka (t a ,x a ,n) dt b B Kb (t b ,x b ,fi)Y\_ dsj J KJ (sj,n) (14) 

\ <^b Wl Qn J 



Here, $Af({gj}) denotes the TV-body massless phase turn is q. The measurement function Mn($n, $l) en- 
space for the N reference jet momenta {qj}, while $l(<z) forces all N jets to be energetic and well enough sepa- 
is the "leptonic" phase space for any additional non- rated so that sy 3> Tn/Q- The index k runs over all 
hadronic particles in the final state, whose total momen- relevant partonic channels, with K a , K b , ■ ■ ■ , kn denoting 
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the individual parton types. 

The hard Wilson coefficient C N is a vector in the ap- 
propriate color space of the external hard partons in each 
partonic channel. It only depends on the hard phase- 
space variables $at and The soft function S N is a 
matrix in the same color space. We can rewrite the color 
contraction as 

C N S N C N = tv(C N C N S N ) = tr(H N S N ) , (15) 

so the hard function Hn = CnC n is also a color-space 
matrix. 

We want to compute the iV-jettiness soft function 



S N (k a ,k b , fci, . . . , k N , {&}, fx) 



(16) 



The ki are the soft contributions to the TJy, so from 
Eq. (fTTj) we have 



P k n s {§3 ■ vk ~ Qi ■ Pk) 



fcEsoft 



(17) 



where the sum now only runs over soft momenta. As indi- 
cated by the second to last argument in Eq. (flo) , the soft 
function still depends on all the reference momenta 
because they enter in the definition of the measured soft 
momentum components in Eq. (|17l) . The soft function is 
defined by the vacuum matrix element 



S$f(k a ,h,ki,...,k N ,{qi}) 



(18) 



= 



Y K H{q i })]j5(k i -2q i .p i )Y K ({q i }) 



where the Pi denotes the momentum operator that picks 
out the total momentum in region i. Here, Y R ({qi}) de- 
notes a product of eikonal Wilson lines in the directions 
in the appropriate path-ordering and color representation 
of the external partons in the partonic channel n. The 
Y K ' and Y K are matrices multiplied in color space. We 
take their overall normalization to be such that the tree 
level result for S N is 1 ]X 6(ki), where 1 is the color iden- 
tity operator (see Eq. (|25p below). 

In the following, we will often use the short-hand no- 
tation 

5^({fci},ju) = S N (k a ,k b , ki,...,kN,{qi},fi) , (19) 

and similarly for other functions that depend on all ki, 
such as the anomalous dimension and counterterm for 
the soft function. 



C. Soft-function RGE 

To derive the structure of the renormalization-group 
equation (RGE) and anomalous dimension of the soft 
function S N (k a , ...,/x), we can use the fact that the fac- 
torized cross section in Eq. (1141) is independent of the 



renormalization scale \x. For this purpose we can ignore 
the phase-space integrals and only consider the last two 
lines in Eq. (fl"4")l . To simplify the notation, we suppress 
the index k and the momentum dependence on the label 
momenta from here on. 

As already mentioned, the hard Wilson coefficient Cn 
is a vector in the color space of the external hard partons, 
so its anomalous dimension jn is a matrix in color space. 
For 1-jettiness (or e + e~ 3-jettiness), the external partons 
are q a (f 'g a , so the only possible color structure is T£p and 
the color space becomes one-dimensional. For q a (f g a g b 
there are already three different color structures = 
{{T a T b ) aP , (T b T a ) a p,8 ab 8 a p}, and so on. 

The hard Wilson coefficient Cn from matching QCD 
onto SCET satisfies the RGE 



M^CV(/j) = 7c(m) C N {p) ■ 
Its anomalous dimension has the general form 



(20) 



7c 



x In 



( _ 1) A„%_^_ i0 l + ^ [asM]) (21) 



where we define A,j = 1 if i and j are both incoming 
or both outgoing partons and A,j = if one of them 
is incoming and the other one outgoing, and with our 
conventions % • qj is always positive. Here, r cusD (q s ) is 
the universal cusp anomalous dimension [30(, 



^cus P (a s ) = -r- 4 - 

47T 



■0(ai 



(22) 



The T" denotes the color charge of the ith external 
particle when coupling to a gluon with color a. It acts 
on the external color space as 



(T?T)... Qi ... = T a Q ift f... ft ... , 

{ L i 1 )-oti- - "Jfta, 1-Pi- 



(23) 



where the first line is for the ith particle being an outgo- 
ing quark or incoming antiquark, the second line for an 
incoming quark or outgoing antiquark, and the third line 
for a gluon. The product 



(24) 



appearing in the first term in Eq. (|2"Tj) . thus represents a 
particular color-space matrix T;j for each choice of i and 
j. We also define the identity operator 1, which acts as 

(if ):::«•;:. = f; (25) 

In particular 

T'j = ld where C q = C q - = C F , C g = C A . (26) 
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With only three partons, q a q^g a , the only possible color 
structure is T 
just numbers, 



structure is T^, so in this case the color matrices are 



1 = 1, = T| = C F 



Tl = C A , 



? • J-g = ~2 °^ ■ 



T • T — T • T 

A 9 ± 9 ~ A 9 A fl 



2 



(27) 



Up to two loops [31| (and maybe more [29j, |32j-|3J|), 
the non-cusp term, 7c (a s ), m Eq. (f2"Tj) is diagonal in 
color and given by a sum over individual quark and gluon 
contributions, 



7c 



(a s ) = 1^7ct(« s ) + O(o?) : 



7c = 7c 



4tt 



3C F + C(^) 



7&(a.) = -gA) + 0(a2). 



(28) 



The RGEs for the beam and jet functions are 

d f 

H—Bi(t,x,fjt) = J dt' -f z B (t - t',fx) Bi(t',x,fi) 



A* ^ Ji ( s ' A 1 ) — / ds ' 7,/( s - s ', A 4 ) J<(s', jti) 



(29) 



The beam and jet anomalous dimension are identical to 



all orders [351 ] . and are given by 

1 / s 

7sO>A*) = 7j(«>m) = -2C i r cusp [a s (/i)] — £ ( — 

A* v A 1 



(30) 



with 



7K« s )=7}(a s ) = -^6C i r + 0(^), 

7> s )=^2/5 o + 0(^), (31) 

and C n (x) denotes the standard plus distribution, 

~0(x)hi n x~ 



C n (x) 



(32) 



Taking the derivative of Eq. (|14p with respect to n, we 
now require 



= [i 



dfi. 



where we use J~i(si, /i) to denote either beam or jet func- 
tions (with s Qi b = t a ,b)i since their RGEs are identical, 
and as before i = a, b, 1, . . . , N. Using Eqs. ([20]) and (|2"9"]) 
together with Eq. (13"3"1) . we get 



+[n*i-».)][*ws»({|}.»')+&({§},M)7 C w + ^&({|}, ( , 



(34) 



where we divided out the Wilson coefficients and shifted the integration variables s$ — > QiT^ — Sj and — > QiTJ| — s[. 
We can now multiply by JJ i J~ 1 {Qiki — QiT^,n) and integrate over TJ^, which replaces JiiQiT^ — s-, /i) — » S(Qiki — 
S0/Qj. Renaming fc^ = Sj/Qi, we obtain 



A^iv({fcJ, = / [II dfc *] ^ " A*) &v({fe-}, M) + aO 7l({fc, - fc*}, A*) 



(35) 



where the soft anomalous dimension is given by 

= 2r cusp [a s ( At )] IE T 4 2 %o (^#) II + E T * • T ^ ln 
- { 1 E 7SK (A*)] + 2 7c [a. (M)] } II S ( k < ) 



c-i^^-ioinw 



= -2r cusp [a s (Ai)]E T * "Tj 



7 ^£o(^S^)+^A ij <5(fc i ) YlSik^ + jsla^YlSih), (36) 

' Sij /X V Sij 
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with the non-cusp part 

7s(a.) = -l£7j(a.)-27&(a.) 

i 

= -lX)[7S(«-) + 27c(«.)]+0(a2) 

= + O(a*). (37) 

In the last step above we rescaled the plus distribution 
X£q(\x) — Cq(x) + In A 5 (x), and applied color identities 
like 



E 



XiTi-Tj = -V(K i +s i )TcTj, (38) 



which follows from color conservation, J7 Tj = 0. 

This derivation shows that factorization implies that 
the kinematic dependence of 7s({fci}, /i) on ki is sepa- 
rable into individual contributions to all orders. This 
generalizes the same result obtained for the special case 
of the hemisphere (i.e. e + e~ 2-jettiness) soft function 
in Ref. (36|, which is reproduced by Eq. (|36|) using 
T t ■ Tj = -C F and 2q t ■ qj = 1. 

Since in Eq. r cusp (a s ), jj(a s ), and 7c (a s ) are 
all known to two loops, so is 7g({fej},/x). The general 
evolution formula in Eq. (1351) leaves Sn hermitian, which 
from Eq. (fT4|) is the only requirement to obtain a real 
cross section. 



Using Eqs. (l3l)|) and (|3T|) . the NLO counterterm is thus 
given by 



Zs({ki},n) 
= lY[S(ki 

i 

i 



<%s(aO i 

2tt 



E T *- T i 



2e 



2 i 



+ yAy (5(fci 



(41) 



Note that since Sn is color diagonal at tree level, the 
imaginary part of Zg does not contribute in Eq. (|39|) at 
NLO, because it cancels between Zg and zt. Hence, 
from Eq. (|4T|) we expect the UV-divergent parts of the 



one-loop bare soft function, S N 3 ' Icl ' 1 \ to have the form 



S 



N 



e(l) 



({**}) = - 



a a (n) 1 



5> 



l+2e 



A- 



n 



+ O{e ) 



(42) 



This implies that the UV divergences are given for any 
N by a simple sum over individual hemisphere contribu- 
tions. We will see how this happens explicitly in the next 
section. 



III. NLO CALCULATION 



D. Renormalization and One-Loop Divergences 



A. General Setup 



The result for the soft anomalous dimension allows us 
to infer the one-loop counterterm for the soft function in 
MS, which we will need in our calculation to renormal- 
ize the bare soft function. This will provide us with a 
nontrivial cross check on our calculation. 

The structure of the anomalous dimension implies that 
the bare and MS renormalized soft functions are related 
by 

S h N arc ({k t }) = J [[[dkl dk'l] ZsiiKh M) (39) 

i 

xS N {{k i -k' i -k'l}, t i)zl{{k'i},ii). 

The bare soft function is independent of /i, so differen- 
tiating both sides with respect to [i determines the soft 
anomalous dimension in terms of the counterterm, 



x (Zs)- i ({k l -k' i }^)^—z s ak' i },^ 



= -2 M A^( {fct '} )M ) + 0( a 2). 



(40) 



Our calculation in the following applies to both 
hadronic and e + e~ collisions, i.e. it is independent of 
whether the Wilson lines are in- or outgoing. For sim- 
plicity we use the term "jet" to refer to both beam jets 
and final-state jets. 

The one-loop diagrams for the soft function are shown 
in Fig. [3l where i and j label any two Wilson lines, and 
we work in Feynman gauge. The virtual diagrams in 
Fig. 3(a) are scaleless and thus vanish in pure dimensional 
regularization. The real emission diagrams in Fig. |3(b)| 
with the gluon attaching to the same Wilson line vanish, 
as it is proportional to qf = 0. Hence, at one loop we can 
write the bare soft function as a sum over contributions 
where the intermediate gluon attaches to the ith and jth 
Wilson line as shown in Figs. 3(c) and |3(d)[ 



S 



barc(l) 
N 



d d p 

Wj 3 (qi-p)(qrP) 



1i ■ 1j 



x2Tr5(p 2 )6(p°)F({k l },{2q l -p}) 



(43) 



The key idea in the hemisphere decomposition method 
is to first fix i and j and then analyze the remaining inte- 
gral. The measurement function resulting from Eq. (|17l) , 
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' % 



(a) 



(b) 





FIG. 3: One- loop diagrams for Sn- The vertical line denotes 
the final-state cut. Diagrams (a) and (b) vanish. Diagrams 
(c) and (d) yield Eq. (|43l . 



which measures the contribution of the gluon in the final 
state to k % , is given by 

F({h}, { P 1 }) = ]T 5(k m - p m ) J] 8(h) 0( P l - p m ) , 

where we denote the component of the gluon momentum 
p M along the jet direction as 



P % = 2% • p . 



(45) 



For example, for 1-jettiness (or e + e 3-jettiness) we 
have three independent labels i ^ j ^ m, so 

F({h},{p 1 }) (46) 
= 6(h - p l ) 6(h) 8{k m ) B(p> - p l ) 6(p m - p<) 
+ 5(h) 5(k 3 - pi) 5{k m ) 6(p i - ^) 6{p m - p>) 
+ 6(h) 5(h) S(k m - P m ) 9(f - P m ) B(p> - p m ) . 

The first two terms correspond to the case where the 
gluon emitted from the «th and jth Wilson line ends up 
in the region of jet i or jet j, respectively. In this case, p 
can become collinear with either qi or qj, resulting in a 
double UV-IR divergence. In the last term, the gluon is 
in the remaining jet m. In this case, both p l > p m and 
pj > pm are bounded from below, so there is only a single 
soft IR divergence. The virtual diagrams, which vanish in 
pure dimensional regularization, turn all IR divergences 
into UV divergences. 

To combine the divergences from the different jet re- 
gions, we split the region of jet m into the two hemi- 
spheres p % < pi and p % > pi defined by the directions of 



jets i and j, 

F({h},{f}) 

= 6(pi - p l ) [5(h - P l ) 5(k m ) 9(p m - p l ) 

+ 5(h) 5(k m - p m ) 9(p l - P m )] 5(k 3 ) + (i o j) 
= 5(h-tf)5(k 3 )6(pi -tf)5(k m ) 

+ [6(h) 6(k m - p m ) - 5(h - P l ) 5(k m ) 
x5(h)6(pi -p l )6(p l -p m ) + (i^ 3 ) 
= Fij^ em i({ki}, {p 1 }) + Fji^ cmi ({h}, {p 1 }) 



F ihm ({h}, {p 1 }) + F ji<m ({h}, {p*}) . 



(47) 



In the second step we replaced 6(p m — p l ) = 1 — 9(p % — 
p m ) in the first term to extend the region for jet i to 
the full p l < pP hemisphere, which gives the hemisphere 
measurement function 

F ijMmi ({ki}, {p 1 }) = 6(h - P l ) 6(h) 6(p> - p l ) 6(k m ) . 

(48) 

The contribution for p m < p l < p> , which overlaps with 
the region for jet m, is subtracted in the second term, 
which gives the non-hemisphere measurement function 
for region m, 

= [6(h) 6(k m - p m ) - 5(h - P l ) 5(k m )] 
x 6(h) 6(p 3 - p l ) 6(p l - p rn ) . (49) 

This splitting up of the measurement function is illus- 
trated in Fig. 21 The integrations over the i and j hemi- 
spheres resulting from -F^hcmi and i*ji,hcmi will now con- 
tain all divergences, while the integration over region m 
resulting from Fij yTn and F 3 i^ n will be UV and IR finite, 
as we will see explicitly below. Essentially, the restriction 
of the emitted and measured gluon to stay away from the 
i and j directions, p IJ > p m , cuts off the UV divergence, 
while the subtraction of the overlapping hemisphere con- 
tribution removes the soft divergence: In the soft limit, 
both p % — > and p m — > 0, and the two terms in square 
brackets in Eq. (l49l) cancel each other. 

We will see in Sec. IIII El (see Eq. (ItTol) below) that this 
split up of the measurement function generalizes to any 
N . Hence, we write the renormalized soft function as 

s N ah}, /i) = i n w + e t * • T ^ 4 l} ({^}' 

i i^j 

+ 0(a 2 s ), (50) 

where we split the NLO contribution into a hemisphere 
and a non-hemisphere contribution, with the latter given 
by a sum over the different jet regions m ^= i,j, 

sS ) ({h},n) = s^l emi ({h},^) + ]T s ( ^ m ({h}^). 



(51) 
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^ TP- ■ ■ TP-- . TP TP 

r ij, hemi r ji,hemi r ij,m r ji,m 

FIG. 4: Separation of the measurement function into hemisphere and non-hemisphere measurement functions for 1-jettiness or 
e + e~ 3-jettiness for a gluon emitted from the ith and jth Wilson line. The phase space is divided into i and j hemispheres 
into which the third jet region m is split. 



Here, Sjj,hemi and Sij >m are the contributions correspond- 
ing to Fjj,hcmi and in Eq. (J47J) . 

It is instructive to compare our hemisphere decompo- 
sition with the method used in Ref. [24| to calculate the 
soft function for cone jets. There the authors first spec- 
ify a jet region and then sum over all contributions from 
different gluon attachments for that fixed jet region. In 
the end they sum over the different jet regions. In this 
case there are nontrivial cancellations between the di- 
vergences (and finite terms) arising from the same gluon 
attachment contributing to different jets. In contrast, as 
seen from Eqs. (|5"0)) and (|5"Tj) . in the hemisphere decom- 
position we first specify a gluon attachment i,j and then 



sum over the contributions to the different jet regions m 
from this specific attachment. This allows us to make the 
cancellations explicit and to isolate the UV divergences 
into the hemisphere contributions. In the end we sum 
over all possible attachments. 



B. Hemisphere Contributions 

Using Eq. (|43[) and restricting the measurement func- 
tion to the hemisphere contribution, i^hemi m Eq. (I47[) , 
we obtain 



-rbare(l) 
r/,hemi 



({ki}) 



IE „2 s 



A 4tt J 



d d p 2sij 



it r(i-e) 

a s (p)l (e^y 



(27r) d p l p> 



2ir5(p 2 ) 0(p°) S(ki - p l ) 5{kj) 9{p> - f) 5{k m ) 



W ' W ?!- ) ^+? - p 1 ) 5{k )e{jP - p l ) 5(k n 



{p 1 p 3 Y 



In the second step we used the coordinate decomposition 



-7' 



P » =p> ^ 



J ij±_ 



(52) 



(53) 



to rewrite the phase-space integral as 



(47T) 



(54) 



and in the last step we integrated over p % and p 1 . The result in Eq. (|52"]) has the expected form in Eq. (l42j) and repro- 
duces the correct counterterm and anomalous dimension. Expanding Eq. (|52[) and subtracting the 1/e divergences, 
we obtain the renormalized NLO hemisphere contribution 



^y, hemi ({**}' A*) ~ J~ 



-A 



Sij (l 



Sij fl 



6 



6(h) 



S(kj) S(k m ) 



(55) 



This generalizes the one-loop result for the hemisphere soft function for two back-to-back jets with equal energies 
from Refs. (3(| to general hemispheres defined by two jet directions g$ and cjj. Note that, as expected from 
reparametrization invariance, the dependence on the jet directions only appears through the invariant sy. 
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C. Non-Hemisphere Contributions 

Wc now turn to the non-hemisphere contributions that account for the precise definition of the 1-jettiness observable 
and the fact that the boundaries between the different jet regions are more complicated than simple hemispheres. 
Inserting the second part in Eq. (|T7|> into Eq. we get 



S 



barc(l) 



ij.ni 



27T 1 



1+e 



5(h) 5(k m - P m ) - S(ki - p l ) 5(k m ) 5(k 3 ) 9(p> - f) 6(f ~ P m ) 



(56) 



To perform the integration over p l we use the rescaled variable x — p J /p l , and rewrite p m in terms of p 1 , x, and the 
angle <f> between q m ± and p± in the transverse plane, 



v' 



p'" 
p' 



1/2 



7771 . Sim ~ ( SjmSim \ ' , / 

H x — 21^-^ — x) cos <p=z(x 

Sin V sr„- / 



(57) 



The limit p m < p l thus implies an upper limit on x, which eliminates the UV divergence for x — > oo. In addition, 
since p m scales like p l , there is also no IR divergence in Eq. (|56"|) , because the term in square brackets vanishes in the 
limit p l — > 0. The integral over p l can then be performed without encountering any divergences, 



M / dp 



0(p l ) 
(piji+a 



S(ki) S(k m - p m ) - 8(ki -p l )S(k m )} = 5(h) -C (^) - -Co (-) 6(k m ) - ln[z(x, 0)] 5(h) 5(k m ) . 

(58) 

Note that the /z-dependence cancels between the first two terms. Taking e —> everywhere else, we obtain the NLO 
non-hemisphere contribution 



)-5(^-£o(— ) 

fj, \ fj,/ U V it / 



ln^5(fci)<5(fc m ) 



<5(/c» ) 5 (fem) 



(59) 



The remaining finite phase-space integrals are defined as (rescaling a; = y 2 (a//3) to simplify the integrands) 



dy 



I (a,p) = - J d0 J jV(y~ #)#-! -j/ 2 + 2j/cos^), 
/ 1 (a,/5) = i y d0 ^ ln(l + y 2 - 2ycos<p)e(y - y/fi/Z) 0(l/a- l-y 2 + 2y 



(60) 



In Eq. (|59]l they are evaluated at a = Sj m /sij and 
/3 = Sim/sij. Their numerical evaluation for fixed a > 
and j3 > poses no problem. We were not able to find 
complete analytic expressions. Their analytic simplifica- 
tion to one-dimensional integrals is given in App.^l with 
the final result in Eq. (|A9J) . 



D. Extension to Other Observables 

As we have just seen, we can extract the divergences in 
the soft function by dividing up phase space into hemi- 
spheres corresponding to pairs of Wilson lines. We will 
now generalize this decomposition to general IR-safe ob- 
servables and to more than three regions (in which case 
the qi are in general non-planar) . 



Consider a measurement that specifies a way to split 
up the angular phase space into non-overlapping regions. 
We use the notation 0j(p) = 1 when the momentum p is 
inside region i and &i(p) = otherwise. We require that 
the union of all regions covers all of phase space, and that 
each region contains at most one of the directions qt, i.e. 

E 8 <(P) = 1 - <\-- (61) 

i 

We explicitly allow the possibility that there are regions 
that do not contain any of the qi, in which case there will 
be more than A + 2 regions. 

In general, we can measure a different observable in 
each region. At NLO, we only need to know how the 
observable for each region i acts on a one-particle state 
with momentum p, which we denote by fi(p). We want fi 
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to be IR safe, which implies that fi(p —> 0) is equivalent 
to measuring no gluon at all. Without loss of general- 
ity we can assume that /i(0) = 0. We will continue to 
denote the arguments of the soft function by ki, which 
are now given by the soft contribution to fi. With this 
notation, the generalization of the measurement function 
in Eq. (|44[) acting on a soft gluon with momentum p is 

F({h},p) = ]T S[k m - f m (p)} e m (p) J] *(*,) . (62) 

We now want to generalize Eq. (|4"T|) by splitting up 
Eq. (|62p into hemisphere and non-hemisphere contribu- 
tions according to which Wilson lines the gluon attaches 
to. We continue to use the labels i and j for the di- 
rections of these two Wilson lines. The hemispheres are 
still determined via the gluon momentum components 
p JJ = 2qij ■ p by p? > p % and p l > pi. Writing out 
Eq. ([52) we now have 

F({h},p) = B(jP - f) [S[h - Mp)} J] 8{h)Q t {p) 

+ S \ k ™ - fmip)] II W e ™W] + (* ° 3) ■ 

(63) 

Note that region j is allowed to overlap with hemisphere 
i, and vice versa. Using Eq. (|6"Tj) . we have 

9i(p) = 1 - £ e ro (p) , (64) 

which allows us to replace the regions i and j by full hemi- 
spheres analogous to Eq. (j47|) . where the complement of 
Qi (p) is now split up between the remaining m (p) with 
m 7^ i. Then Eq. (|52^) can be written as 

F({h},P) = F ljMmi ({ki},p) + Fj iMmi ({ki},p) (65) 

+ X! F i3,m{{h},P) + Y l\,,.,n ■ {>'■,) ■ I" , 

where the hemisphere contributions are given by 
F ijMmi ({ki},p) = 6{p? -p*) 5[h - fi(p)} ]j5(ki) , (66) 

and the non-hemisphere contributions by 

Fij, m ({ki},p) 

= 9&-p i )Q m (p) TJ 5(h) (67) 

x {6 (hi) S[k m - f m (p)} - 5[h - fi(p)\ S(k m )} . 

As in Sec. IIII C[ all the divergences are contained in 
the hemisphere contributions, while the non-hemisphere 
contributions are UV and IR finite. The measurement of 
either fi or f m in Eq. (|67j) fixes the magnitude of p, while 
the restriction of the emitted gluon to region m forces 



it to stay away from the i and j directions. Taken to- 
gether this eliminates the UV divergence. The IR-safety 
of fi then ensures that in the limit p — > the terms 
in curly brackets in Eq. (|67[) cancel each other, which 
eliminates the IR divergence. As a result, for any set of 
IR-safe observables fi all UV divergences, and hence the 
anomalous dimension, are contained in the hemisphere 
contributions determined by Eq. (|66p . Depending on the 
observable, these contributions can be more complicated 
than in Eq. (I5"5j) . Note that this result depends on the 
fact that an observable fi is measured in each region i. 
If we have a region u where only an angular restriction is 
imposed by & u (p), the corresponding S[k u — f u {p)) is ab- 
sent (an "unmeasured jet" in the language of Ref. ^23]). 
In this case the hemisphere contributions S'uj.hemi are 
scaleless and vanish. The non-hemisphere contributions 
S u j,m and Sij tU are still IR-fmite, but now contain a UV 
divergence in the term coming from region u, for which 
the magnitude of p is not fixed anymore. In this case, the 
factorization structure is different and the soft anoma- 
lous dimension depends on the parameters determining 
the boundary of region u, for example the cone radius as 
in Ref. 

Although we have only applied the hemisphere decom- 
position method at NLO, the V-jettiness factorization 
theorem implies that the UV divergences and soft anoma- 
lous dimensions factor into pairwise hemisphere contri- 
butions to all orders, as shown by Eq. (f36|) . Hence, we 
believe the hemisphere decomposition will remain useful 
also at higher orders. 

E. NLO Calculation for AT-Jettiness 

We now use the general arguments in the previous sub- 
section and apply them to the case of V-jettiness. In 
this case the observables are simply the components of 
the gluon momentum along the jet directions, while the 
regions are determined by the smallest p % . Hence, 

f t ( P ) = p i = 2q i -p, Qi(p) = TJ 6(p m ~ p l ) , (68) 

which turns Eq. (1521 into Eq. flUJ. From Eqs. (|6"6")l and 
([ST) we get 

F ijMai {{ki}, {p 1 }) = 9(p> - f) S(h - f) Y[ 5{k m ) , 

(69) 

and 

F lhm {{h},{p 1 )) (70) 
= [5(h) S(k m - p m ) - 5(h - p l ) S(k m )] 

x 9(p> - p l ) 9(p l - p m ) Yl 6(h) 6{p l - p m ) . 

The calculation of the hemisphere contribution for gen- 
eral V is identical to the 1-jettiness case in Sec. IIII Bl with 
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the overall replacement 6(kj) S(k m ) —> Y\ m ^i5{k m ) aris- 
ing from Eq. (|69j) . In particular, we can see immediately 
that this reproduces the correct NLO counterterm and 
soft anomalous dimension in Eqs. (|4T|) and (|36|) . The fi- 
nal result for the renormalized hemisphere contribution 
is given by Eq. 



a a (n) 



47T 

ze! 

6 



-A 



6{k 



n <km 



(71) 



For the non-hemisphere contribution, there are now 
several regions m contributing. The calculation for each 
region proceeds as in Sec. lIII (Jl except that we now have 



additional 8(p l —p m ) functions in Eq. (|70|) . which separate 
region m from the remaining regions I ^ i,m. We can 
write p m and p l in terms of p l and x = pP /p l , 



p m 


_ Sj m 


. Sim 


p' 


Sij 


Sij 


t 


= f£ 4 


Sil 








p' 


Sij 


Sij 



.fSjmSim \ , 
1^^-72 Xj COS0 



1/2 



(72) 



Here <fi is again defined as the angle between p± and q m ± , 
while (f>i m are the angles between the remaining qi± and 
q m ±. The result for S^ m ({ki}, ^) has the same form as 



Eq. ([591), 



7T I \ Sin Sii LS™. S im U^i,j,m/ U V /i / U \ U 



In -r — 



5(fci) 8{k m ) 



[1 V/i. 



llife^i^^.d.. )<W(u) n TO- (73) 

V Sii Sii *~ Sim Sim. ' I _ 



l^i.m 



The finite phase-space integrals are now given by 

'dy 



I (a,/3,{ai,Pi,4>i}) 



(y- ^JJ^)6(l/ a -l-y 2 + 2yt 



x J] e\ai - 1 + {fa - l)t/ 2 - 2y[^fa cos(0 + 00 
2 



ii(a,/3,{a z , = - 



— ln(l + y 2 - 2ycos0) #(y- ^/a) 0(l/a- 1 - y 2 + 2y cos< 



x Y[6[ai - 1 + (fa - l)y 2 - 22/fv/a^ft cos(0 + </>;) - cos< 



(74) 



An algorithm to systematically evaluate them numerically is given in App. [X] The values for the parameters in 
Eq. (O are a = Sjm/kj, P = Sim/%> "2 = Sjl/sjm, and fa = su/s im . 



IV. CONCLUSIONS 

./V-jettiness is a global event shape that can be used to 
define an exclusive iV-jet cross section. We have given a 
factorization theorem for the cross section fully differen- 
tial in the individual ./V-jettiness contributions for each 
region, TJ^, which correspond to the mass of each jet re- 
gion. We have computed the corresponding ./V-jettiness 
soft function, differential in all TJ^ at one loop. 

In our calculation we analytically extract the UV di- 
vergences by splitting the phase space into hemispheres 
depending on which Wilson lines the soft gluon attaches 
to. The hemisphere contributions reproduce the anoma- 
lous dimension of the soft function as expected from the 
consistency of the factorization theorem. The remaining 
non-hemisphere contributions, which encode the depen- 



dence on the boundaries between the regions, are reduced 
to one-dimensional numerical integrals. We show that 
this hemisphere decomposition can be applied in general 
to compute soft functions for other observables, such as 
jet algorithms and jet shapes, at one loop. We also expect 
that it can be generalized to two loops. 

Our soft-function calculation provides the last missing 
ingredient to obtain the exclusive -/V-jet cross section re- 
summcd to NNLL for any process where the correspond- 
ing SCET hard function at NLO is known from the one- 
loop QCD calculation. In many processes it has been ob- 
tained explicitly [28l [38l445j . In general, the NLO hard 
function is given in terms of the virtual one-loop QCD 
diagrams, and there are large ongo ing, efforts to compute 
these for many LHC processes pa - l54| . 

The shape of the jet regions as determined by N- 
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jettiness depend on the specific distance measure used, 
and our results apply to any choice of distance measure. 
As we saw in Fig. [2j using a geometric measure, the jet 
regions yield jets with circular boundaries, which is a fea- 
ture desired experimentally. Hence, it will be interesting 
to explore the use of iV-jettiness directly as an exclusive 
iV-jet algorithm in the future. 
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Appendix A: Finite Integrals 

1. 1-Jettiness 

Here we further study the finite phase-space integrals 
in Eq. (|60]l that are required for 1-jettiness or e + e~ 3- 
jettiness. The indefinite integrals over y can be carried 
out explicitly. In particular, for 1\ we have 

G(y,tt>) = J^Hl+y 2 -2ycos(ft) =-2Re[Li 2 (ye i *)] . 

(Al) 

The remaining integrals over (ft must be be carried out 
numerically. 2 

What remains is to determine the ^-dependent inte- 
gration limits on y. We use — tt < (ft < ir as the funda- 
mental region for (ft. Also recall that a = §j m /sij and 
P = Si m /sij, which are positive definite. The 8 functions 
in Eq. (IBTfl) impose the conditions 



(y — cos ( 



sin 



1 

< - . 

a 



(A2) 



which are illustrated in Fig. [5] Solving for y they imply 



, a) f < V < V+ [<ft, ") , sin (ft < - , 

a 




y_ (cf), a) = cos (ft — \j\ja~ sin 2 , 
y + ((ft,a) — cos (ft + yl/a — sin 2 . (A3) 
We can now distinguish the two cases a < 1 and a > 1. 



2 One could also think about first integrating over <f>, since the 
original ^-integral can be done and the limits are linear in cos <j>. 
This does not lead to any simplification, however, because the 
remaining numerical y-integral will then involve arccos[(l + y 2 — 
V«)/(2»)]. 



Case a < 1 For a < 1, we have Sj m < s^, which 
means that cjj is closer to q m than to cji. In this case, 
which is illustrated in the left panel of Fig. O the roots 
always exist and y~(<p) is strictly negative, so we have 



< y < y+0,a) 



-a + (3 - 1 



1 > cos</) > max|^— r — == — . — 1 



< Va + 1 



(A4) 



The lower limit on cos <f> is necessary to guarantee that 
y/P/a < y+ (4>, a) ■ The condition on a and (3 is then 
necessary to guarantee that 1 > (a+/3—l) / (2y/af3), such 
that the lower cos <j> limit does not exceed the upper one, 
otherwise the integral vanishes. The lower cos</> limit 
itself is only nontrivial for <Ja + \fp > 1 which means 
\Jsim + J Sj m > \Jlkj. For a purely geometric measure 
this is always true, but it need not be the case for more 
general measures. 

Case a > 1 For a > 1, illustrated in the right panel 
of Fig. [SJ the condition sin 2 (j) < 1/a for the roots to exist 
becomes nontrivial and forces an upper limit on \<j)\, 



\<fi\ < arcsin —= . 



(A5) 



(The second solution for the arcsin is not allowed for 
y > 0.) Now, both lower limits on y are possible. To 
determine which y-limit applies at a given value of (ft, we 
can distinguish two cases. First, 



- < y < y+(4>,a) , 

a 



1 > COS (ft > 



q + jg- 1 
Vu- 1 < \fii< Va+1, 



(A6) 



where the cos (ft limits result from enforcing y_ < 
y/P/a < y+ and the conditions on a and P enforce the 
lower limit on cos<f> to be < 1. Second, 



y-{<t>,Oi) < y < y+{(ft,a) 
a + P-1 

2y/^P 



> COS (ft > 



1 

1 — . 

a 



P<a-1 



(A7) 



where the upper cos (ft limit and the condition on a and P 
arises from requiring y/P/a < y^((ft,a), while the lower 
limit on cos (ft is equivalent to Eq. (|A5[) . 

Combined Result To write the various conditions in a 
compact form we define the following two angles 



0ma X (a) = arcsin 



1 



la 
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-e- o 



I 

(y — cos 


a < 1 
4>) 2 + sin 2 4> < 1/a 







^cut 

-e- o 







a > 1 



(y — cos'(/>) 2 + sin 2 </> < 1/ot 



FIG. 5: Phase-space constraints from Eq. (|A3[l in the (\>-y plane for a — 0.8 (left) and a — 1.2 (right). In both cases /3/a — 0.5. 



4>cut{a,(3) 



arccos ■ 



a + 0-1 



> i, 

v^ + V?<i, 

otherwise . 

(A8) 



The conditions for a < 1 and the first case for a > 1 
reduce to |0| < </> C ut- For the second case for a > 1, which 
only applies for /? < a — 1, we have cu t < \<f>\ < </> m ax- 
Using the fact that the integrand is symmetric in </>, the 
final result for the integrals is given by 



I {a,P) = 2 — In — -==- + W{a - fi - 1) / — In — — — -, 

Jo n \/P/a J0 cut (Q,/3) T y-{9,a) 



h(a,0)=2 / ^[G(y + ^,a)i)-G(vWa^)] 
Jo ^ 

+ 20(a-/3-l) r°" X( ^ [G(y + (0,a),0) -G(2/_(0,a),0)] . (A9) 
2. JV-Jettiness 

We now turn to the integrals Io,i(o<, /3, {aj, /3j, c/f>z}), defined in Eq. (|74")l . that are needed for general TV. The y- 
intcgral is the same as before and can be carried out explicitly. For a given value of </>, the 9 functions split the y 
integration region into a number of mutually exclusive y-intervals, which yields 



J (a, 13, {a u j3 h <M) = f ^ £ln %^) - ^ in (0)] , 

h{aMa u Pu<k})= I — E[ G ^n 1 a X (0)^)"G(yi i j0),0)]^[yL x W)-yii„(^)]- (A10) 



7T 



Here, the sum runs over all intervals and 2/4in(<W an d Um&xi't') are the lower and upper limits of the Jth interval, and 
can depend on all a, /3,ai, j3i,<pi. 

What remains is to determine the y-limits for a given (f>. The conditions imposed by the primary 9 functions 
involving a and /? are as in the previous subsection. The additional 9 functions impose the condition for each I 



1-ai + (1 - Pi)y 2 - 2y[cos0- ^/arfh cos(0 + 0j)] < °- 



(All) 
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Recall that ai = Sji/sj m > and pi = su/si m > 0. They essentially compare the distance between qi and qtj with 
the distance between q m and qi.j. The angle <pi = 4>lm is the angle between qi± and q m ±- The limits on y coming 
from Eq. (| Al 1|) are given in terms of the roots of the polynomial, 



y±(<£,aj,/3j,(fo) = ^ cos0- cos(0 + z ) ± y [cos0- a/ ft cos(?!> + 0;)] - (1 - a ; )(l - A) > . (A12) 



To analyze the limits on y imposed by Eq. (IA11[) for 
each I, there are three questions to ask: 

1. Does the parabola open upwards or downwards? 

2. Does it have real roots? 

3. What are the signs of the roots? 
The condition for the roots to exist is 

[cos 4> - V^cos(0 + <h)] 2 > (l-ai)(l-A)- (A13) 

The correct y limits at a given fixed value of <f> are then 
determined as follows: 

1. Pi < 1: The parabola opens upwards, so y must be 
in between the two roots, y_ < y < y+. 

(a) oli > 1: Equation (IA13[) is always satisfied, 
V- < 0, and y+ > gives an upper limit 

y<y + ((f>,ai,l3 h (f>i). (A14) 

(b) ai < 1: Equation (|A13|) is nontrivial, and the 
roots have the same sign if they exist. Hence, 

y-(<P,oti,0i,<pi) < y < y+((j>,ai,Pi,<l>i) , 
cos</>- y/atfhcasty + fa) > (A15) 

The y-integral vanishes if the condition on <j) 
is not satisfied. 

2. /3i > 1: The parabola opens downwards, so y must 
be outside the two roots, y < y_ or y > y+. 

(a) aj < 1: Equation (IA13[) is always satisfied, 
V- < 0, and y+ > gives lower limit 

y>y+(<t>,ai,Pl,<t>l)- (A16) 

(b) oli > 1: Equation (|A13[) is nontrivial, and the 
roots have the same sign if they exist. Hence, 

V <y-(<P,ai,Pi,(f>i) or y>y + {cj),ai,Pi,(j>i), 
cos^- v/q^Acos^ + <£,) > v/(i-az)(i-A)- (A17) 

There arc no constraints on y if the condition 
on cj) is not satisfied. 



I 

3. Pi = 1: There is no parabola. 

(a) a; < 1: The limits are 

. 1 - oti 

y — 1 

~~ 2cos0 — 2^/ai cos(cj) + 0;) ' 
cos 4> > v^ol cos(</) + </>/), (A18) 

and the y-integral vanishes if the condition on 
4> is not satisfied. 

(b) ai > 1: The limits are 

. ai ~ 1 

2y/ai cos(0 + 0; ) — 2 cos </> 
cos 4> < v^cty cos((/) 4- 4>i) . (A19) 

There arc no constraints on y if the condition 
on <j> is not satisfied. 



In principle one can now combine all limits and de- 
termine all possible ^-intervals in which a particular 
set of lower and upper y-limits applies, as we did in 
Eq. (|A9[) . However, although this is straightforward it 
quickly becomes very cumbersome. Alternatively, it is 
easy to devise an algorithm to obtain the correct y-limits 
m Eq. ([ATO]) for a given value of <f> in the numerical inte- 
gration over 4>. One starts with the y-limits in Eq. (|A3I) . 
call them y m - m and y max - Next, one loops over all I and 
determines the limits imposed by each I as above. If one 
encounters a stronger lower or upper limit, y mm and/or 
2/max are updated to the new stronger limit. If one en- 
counters a necessary condition on <fi that is violated, the 
integrand vanishes and one can stop. Case 2(b) requires 
special attention. If it is encountered, the y interval is 
split in two if necessary and one continues by maintaining 
two (or more) mutually exclusive y-intervals each hav- 
ing its own lower and upper limit. Newly encountered 
stronger limits are then applied to each interval. An in- 
terval is eliminated whenever its lower limit exceeds its 
upper limit. If the last existing interval is eliminated the 
integrand vanishes. 
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